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We present a novel derivation of the Teukolsky-Starobinsky identities, based on properties of the 
confluent Heun functions. These functions define analytically all exact solutions to the Teukolsky 
master equation, as well as to the Regge- Wheeler and Zerilli ones. The class of solutions, subject to 
Teukolsky-Starobinsky type of identities is studied. Our generalization of the Teukolsky-Starobinsky 
identities is valid for the already studied linear perturbations to the Kerr and Schwarzschild metrics, 
as well as for large new classes of of such perturbations which are explicitly described in the present 
article. Symmetry of parameters of confluent Heun's functions is shown to stay behind the behavior 
of the known solutions under the change of the sign of their spin weights. A new efficient recurrent 
method for calculation of Starobinsky's constant is described. 

PACS numbers: 04.70. Bw, 04.30.-w., 04.30 Nx 



I. INTRODUCTION 

Today the Teukolsky-Starobinsky identities (TSI) are 
essential ingredient of the theory of perturbations to a 
gravitational field of Kerr black holes. For spin one and 
two the TSI were discovered together with their first ap- 
plications in the study of perturbations of rotating rela- 
tivistic objects pioneered by Teukolsky [l[ and Starobin- 
sky • An independent derivation and detailed study of 
TSI was given by Chandrasekhar in [U, Q and by other 
authors in |5j. Later on a proper generalization of the 
TSI for all physically interesting spins, as well as the TSI 
in the presence of a nonzero cosmological constant was 
found and justified in [f|. Some more recent applications 
can be found in 0]. 

All known derivations of TSI are based on a direct 
use of the famous Teukolsky master equation (TME), or 
its generalizations, and on some very special properties 
of the solutions in use. For example, one utilizes the 
regularity and integrability of the used angular solutions, 
parameters with real values, etc. 

The TME describes the perturbations ty(t,r,6,tp) 
of fields of all interesting spin weights s — 
0, ±1/2, ±1, ±3/2, ±2 in Kerr background in terms of 
Newman-Penrose scalars. Various significant results and 
references may be found in [l|, 0, H, 0, H, H, 0, EH • 

The key feature is that in Boyer-Lindquist coordi- 
nates one can separate the variables using the ansatz 
ty(t,r,8,(p) = exp(—iujt + ini(p)R(r)S(9). Thus, a pair of 
two connected differential equations arises - the Teukol- 
sky angular equation (TAE): 
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and the Teukolsky radial equation (TRE): 



A" s |- (V^i?^)) + (V(r) - fi)iZ(r) = 0, (1.2) 



V(r) = Msior + 2mato 
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The real parameter a = J/M > is related with angu- 
lar momentum J of the Kerr metric (a < M for black 
holes, or a > M - for naked singularities), M is the Ke- 
plerian mass of the Kerr solution, K — uj(r 2 + a 2 ) — ma, 
A = (r — r + ) (r — r_ ) = r 2 — 2Mr + a 2 . The two complex 
separation constants lu and E are to be determined using 
the boundary conditions of the problem. 

At present the usage of the Kerr metric for the de- 
scription of astrophysical black holes is widely accepted. 
The astrophysical application of the naked singularities 
seems to be more problematic because of their instability 
P. However, such instability may be useful for construc- 
tion of models of astrophysical explosions like gamma ray 
bursts Q . Our method for the derivation of the general- 
ized TSI is applicable to both cases, despite of the fact 
that some quantities (like r±) are complex for the naked 
singularities. 

The TSI for angular Teukolsky's function S S(9) of ar- 
bitrary spin weight s can be written in the form 

2 _ s iL S L( +S S(9)) = s £(- s S(e)), (1.3) 

i-.tfa-.tf ...^tf s tf (sS(9)) = S €( +S S(9)); 
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Using the operators 
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one obtains the TSI for the radial Teukolsky function 
s R{r) of arbitrary spin weight s in the form [181 ]: 



R(r)) 

A* ( Dt)" (_ s i?(r)) 



A s ( D)" S (A s +S 

2s 



.2)(_,iJ(r)), (1.6) 
S D*(A S +s R(r)). 



In Eqs. p.3p and (JT76J) ^ and s 2) are the corresponding 
Starobinsky constants The standard calculation of 
these constants for arbitrary s is quite complicated [!, H| . 

In the widespread derivation of TSI, proposed for the 
first time in Q , one supposes to work with real frequen- 
cies uj. This is explicitly stressed in 0], as well as in some 
of the articles [1,0. As a result, the operators T) and 
can be considered conjugated; this property is system- 
atically used for derivation of the TSI d, 0, HI, @ • How- 
ever, the well-known boundary problems, related to black 
holes and naked singularities, yield complex frequencies 
lu with certainly nonzero imaginary part 3(o>) ^ (Se e 
for example the Refs. [3, [l(| and the recent articles [l(| , 
where a new boundary problem was studied and some 
very preliminary results of our attempts to apply singu- 
lar solutions of TAE to the description of collimated rel- 
ativistic jets are reported.) The physical reason is clear: 
The systems under consideration are open physical sys- 
tems and their energy is not conserved. Hence, for them 
Q(cj) 7^ 0. Our derivation of TSI does not imply complex 
conjugation of the corresponding quantities. It is valid 
for complex frequencies u>, as well as for complex values 
of all other parameters. Despite the fact that our deriva- 
tion is quite different in form, it is closely related to the 
original one, outlined in [H, 0|. This original method is 
also valid for complex frequencies. 

It has been well known for a long time [ll| that for 
a^M 20] the TAE (H) and TRE fL2]) can be reduced 
to the confluent Heun ordinary differential equation [l5| . 
Recently in [l2j all classes of the exact solutions to TAE 

1|) and TRE (|I.2[) were described using the confluent 
Heun function HeunC(a, 0, 7, S, rj, z) - a unique particu- 
lar solution of the confluent Heun equation, written here 
in the simplest uniform shape: 



z-l 



z z — 1 



H = 0. (1.7) 



The function HeunC(a, 0, 7, S, rj, z) is defined as a solu- 
tion of Eq. (|I.7p which is regular in vicinity of the sin- 
gular point z = and subject to the normalization con- 
dition HeunC(a, 0, 7, 5, rj, 0) = 1 [Uj]. The parameters 
a, 0, 7, 6, rj, introduced in [151 ] and used in the widespread 
computer package MAPLE, are related with /1 and v, ac- 
cording to the equations /j,— ^(a~ — 7 + a0 — 0j) — rj 
and v = i(a + + 7 + a-f + 0-f) + S + rj. 

The other particular solutions to Eq. (|I.7[) are not 
termed "confluent Heun's functions" , according to the ac- 
cepted modern terminology [lj| . The reason is that in the 
general case other solutions can be represented in a non- 
trivial way in terms of solutions HeunC(a, /?, 7, 5, 77, z). 
An example if the application of this specific property of 



the Eq. (|I.7[) is the basic formula pi.l[) in the next sec- 
tion. Therefore from a computational point of view it is 
sufficient to study the Taylor series of this standard local 
solution and its analytical continuation in the complex 
plane C z . Thus, the instrumental use of the confluent 
Heun function HeunC(a, (3, 7, 5, 77, z) is much more ad- 
vantageous than the simple fact, recognized already in 
[iH ]. that the TRE and TAE can be reduced to the Eq. 

Essential novel properties of the confluent Heun func- 
tions HeunC(a, (3, 7, 5, 77, z) were studied for the first time 
in the recent article [13[. These properties are prerequi- 
sites for our derivation of the TSI in the next section. 
The subclass of <5Ar-confluent-Heun's-functions 

HeunCAr(a, 0, 7, 77, z) = HeunC(a, (3,j,8n, 77, z) , (1.8) 

was introduced assuming the fulfilment of the "<5at- 
condition," 



fo = -a(JV+l+09+7)/2), 

for some fixed nonnegative integer N = 0, 1, 2, . . . 
For functions (11.811 and their associate functions 



(1.9) 



HeunC* (a, /?, 7, 77, z) = (1. 10) 

= HeunC (a^ , , 7* , , 77* , z) , 



a*=a, I3% = (3+N+1, 7* = 7 +iV+l, 



(1.11) 



5J L= 5_ 
a% a 



+N + l, V % = ri- 



(iV+l)(AT+l-x) 



x=a—[3-^i] the following simple basic relation is available: 
d N+1 

dz N+i HeunC w (a, /?, 7, 77, z) = (1.12) 
= *Pjv HeunC* (a, 0, 7, 77, z). 

¥ N = {N+l)\v N+1 (a,l3^ 1 -a{N + l + ([3+ 1 )/2),r ] ) is a 
constant related to the coefficients in the Taylor series 

00 

HeunC(a,/3,7, S, 77, z) = ^ v n{&, f3, 7, S, 77)2". (1.13) 

n=0 

These new results are the mathematical basis of the 
present article. We derive TSI using the above properties 
of the confluent Heun function HeunC(a, 0, 7, 8, rj, z) and 
show that the Starobinsky constant is related to its Tay- 
lor series coefficients. This yields a new effective method 
for the calculation of Starobinsky's constant for any spin. 

Our derivation of the TSI, being simpler, universal and 
straightforward, supplements and justifies the usual ap- 
proach. It is uniform and valid for both the TAE and 
TRE. It extends the TSI to new classes of solutions to 
the TME, described below. Thus the known fascinating 
properties of the standard solutions of the TME acquire 
a natural mathematical framework and proper extension. 
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Moreover, the same derivation yields Teukolsky- 
Starobinsky like identities for the Regge- Wheeler equa- 
tion (RWE) and Zerilli equation (ZE). To our knowledge, 
until now the TSI were not known for the RWE and ZE. 

Our approach is in accord with the Chandrasekhar ex- 
pectation that the Newman-Penrose formalism for per- 
turbations of the Kerr metric may "...enable us to dis- 
cover new classes of identities among the special func- 
tions of mathematical physics when they occur as solu- 
tions of Einstein's equations" and his suggestion to " ...in- 
clude Teukolsky's functions among " special" functions of 
mathematical physics..." 3]. Today it is clear that for 
realization of this general idea we need to use the Heun 
functions and some of their generalizations [l7j |. 



II. THE EXACT SOLUTIONS TO THE TAE, 

TRE AND RWE IN TERMS OF THE 
CONFLUENT HEUN FUNCTION AND THE 
GENERALIZED TSI 

The exact local solutions around the singular points 
z± = z T = 1) of the three differential equations 
listed in the introduction - the TRE, TAE, and RWE 
- can be written in the following common and universal 
form pi]: 

s^ E , m , {<Ty (z ± )= (HI) 

= g(z ± ) e az ±/ 2 z m z 7/2 HeunC(a, 0, 7, 5, r), z ± ). 

The TSI turn out to be a representation of the basic re- 
lation (|I.12|) in terms of the corresponding (^-solutions of 
type (jll.ip . Indeed, let us use the short notation X^(z±), 
X*(z±) for the <5/v-solutions (|II.1[) which correspond to 
the functions p.8[) . (|I.10[) of argument z± in the place of 
the factor HeunC(a, (3, 7, 8, 77, z ± ) in Eq. (|II.1[> . Then, 
using relations p.ll|) and g* — g, we obtain from (|I.12|) 
the following universal form of the generalized Teukolsky- 
Starobinsky identities: 



( 



: z_) {N+1)/2 D% +1 X N (z ± ) =y N X*(z ± ), (11.2) 



where 



^ = ^-?t|^±J^--^. (H.3) 



dz± 



z± 



g dz± 



For compactness, in the right-hand-side of Eqs. 1|) . 
(III.3P we are using the following abbreviated notation: 
a = a a a±, (3 = (Jf 3 f3±, 7 = ct 7 7±; {a} = {0^,073,0^}, 
&a,f3,~t = ±1- Thus, we obtain 16 local solutions to any 
of the TRE, TAE and RWE, making use of the confluent 
Heun's function HeunC(a, /3, 7, 5, 77, z) of different argu- 
ment z+ or z_ in the right-hand-side of Eq. (jll.lj) and 
changing the signs of the parameters a,/3, 7. 

The difference among the solutions to the above three 
problems (the TRE, TAE and RWE) is in the values 
of the parameters and arguments in their common form 
(ULU: 



1. According to [II], for the TAE fllT) q(z ± ) = 1 and 
a± — ±4aw, (3± —s =F to, j± =s ± to, 5± = ±4saw, 
v± = !^!±f! T 2s auj -aW-E; (11 A) 



z+ = z + (6) = (cos(6>/2)) 2 , z_ = Z-{6) = (sin(6»/2)) 2 . 

2. For the TRE (|L2]) g(z ± ) = A" s / 2 . Using the quan- 
tities p= y/r+/r--y/r-/r + , tt± = y/r T /r±/ (r + +r_), 
n a = (n + + Q-) /2 = l/(2a), n g = ^/n + n_ = l/(2M) we 
have: 



a ± — ±2iuj(r + — r_) = ±ipuj/Q, a , 

P ± = S ± 2i (TO — Uj/n^:) /p, 

7± = s T 2i (to - /p, 



a± (s - icu/flg) , 



(II.5a) 
(II.5b) 
(H.5c) 
(II.5d) 



2 1 2 

S + TO - 



-E 
(2mfl a - lo) 2 



2m 2 n 



202 



p^ni 



1 



2 V ^a^g t 



-, Z- ■ 



(II.5f) 



3. For the Schwarzschild metric the discussed results 
(in terms of the Weyl scalars) can be obtained from the 
case of the Kerr metric in the limit a — > 0, see [Hj]. 
An independent treatment is possible, ma king use of the 
RWE. Then in units 2M = 1 we have [HEif: 



g = r, E = 1(1 + 1), I = \s\, \s\ + 1, .. . ; (II.6a) 

2iio 
Is 



(II.6b) 



S ± = ±2lu 2 , 7 1± 

z + = r, Z- = 1 — r 



-E + s 2 
-E + s 2 + 2lj 2( ' 



(II.6c) 
(II.6d) 



The knowledge of the 16 local solutions (jll.ip to the 
each of the TRE, TAE and RWE makes it possible to 
write down their local general solutions as linear combi- 
nations of different pairs of linearly independent particu- 
lar solutions. Then we can formulate different boundary 
problems for the TRE, TAE and RWE For example, 
one can impose the standard black hole boundary con- 
ditions on the solutions of the TRE and combine them 
with the standard regularity requirement of the solutions 
to the TAE on the poles 9 = 0, ir. This way we obtain 
the standard quasinormal modes of the Kerr black holes 
in terms of confluent Heun's functions [l2l ]. Considering 
polynomial solutions of the TRE and TAE we are able 
to describe collimated one-way running waves of pertur- 
bations to the Kerr metric in terms of confluent Heun's 
polynomials [l2j]. One may hope that the proper usage 
of the singular solutions to the TAE may help us to de- 
scribe the creation of relativistic jets by such collimated 
running waves, see [T2I [l||. These important physical 
problems need a separate careful investigation. 
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In the present article we are concentrating our efforts 
on the mathematical properties of the solutions related 
with the TSI. Knowledge of these properties is necessary 
for further progress in the field. As we shall see in the 
next sections, it is natural to formulate these common 
general properties in terms of confluent Hcun's functions. 

III. 5jv-CONFLUENT HEUN'S FUNCTIONS 
AND SOLUTIONS TO THE TAE, TRE AND RWE 

The fulfilment of the Sn condition for solutions to the 
TAE, TRE and RWE is tigh tl y rel ated to th e specific 
values of the parameters (|II.4|1 . (|H.5|) . or (|II.6[) . One can 
find a detailed general consideration of this problem in 
0. 

For the TAE the Sn condition yields two different 
cases: 



which the J^-condition yields important extra restric- 
tions on the free parameters of the problem. We call 
"^-solution" any solution (jll.lj) to the TAE, TRE, or 
RWE, subject to the <5jv condition p.9[) . In section 4 we 
will use the term " (^v-solution" for the solutions (|IV.6b|> 
to the ZE, having in mind the fulfillment of the Sn con- 
dition (|L9|, too. 

Up to now only ^-solutions to the TRE and TAE 
which obey the relations (|III.1[) were studied. The Sn- 
solutions, which obey the other relations - ()III.2[) . for 
the TAE, (|IlT3l) - (|nL4l ). for the TRE and (|IIL6| . for the 
RWE, are new. Their physical meaning and applications 
still have to be recovered. In the next section we derive 
the generalized TSI, valid for all Sn solutions. 



IV. A NOVEL DERIVATION OF THE 
TEUKOLSKY-STAROBINSKY IDENTITIES 



s N+l = 2\s\ > 1, |s| > 1/2, a a = ai3 = a^ = -a.(lll.l) 



s N m + l = \s\ ± ma > 1, a a = — ap = er 7 = — a. (III. 2) 

For the TRE the Sn condition yields four different 
cases. In the first one relations 1 1) are also valid. 
Note that there a — sign(s). In the other three cases 
instead of the integer N, the frequencies s uj n m CTq a a ^ 
are fixed: 



±i^±i (III.3) 
4M ' V ; 



(III.4a) 



(III.4b) 



= mfi_ ± — ( — - 1 ) (N + 1 ± s). 



The S n condition yields the following cases for the 
RWE: 



For all eight classes of <5Ar-solutions to the TRE, TAE 
and RWE, described in Sec. Ill, one can derive in a uni- 
form way identities of the Teukolsky-Starobinsky type, 
starting from the basic form (|IL2|) . Relations (jlL3|) . (jTTT4j) 
and (|II.5|) yield the following explicit representations of 
the generalized TSI [H|: 

1. In the case of the TAE, 



/ - \ N+l 

^ 6 ) N+1 (& s N (e)=£ N s N *(ey 



,± = — + cr a aLU sm( 
at) 



2sin6» 



(^ + (7 7 )sT(^-a 7 )m 
2 

S% (9) = X N (z±), 
Cn = (±1) N+1 ¥n. 



(IV. la) 



(IV.lb) 

(IV. lc) 
(IV.ld) 



Disentangling the differential operators in Eq. pV.lap , 
we can represent it in a form similar to Eq. (|I.3[1 : 



N+l 



11 (d e ,± - kcot&)s%(9) = t N S%*(9). (IV.2) 



s N+l = \s\ > 1, 

<T Q =(7-y, <7p = — (J, Or <7 a =ap, (7-y = —(J, 



(III.5) 



{N+l + aps), a a = -a p ; (III.6a) 



Here an ordered operator product is used. The arrow in- 
dicates the operator ordering and points in the direction 
of the increase of the integer k. 
2. In the case of the TRE, 



^ = l -^(N+l+a J s),<r a = -^. (111.6b) 



Until recently only solutions of a very special kind were 
studied and used in the existing literature. For these spe- 
cial solutions the relations 1|) take place and the Sn~ 
condition is fulfilled automatically and simultaneously for 
both the TRE and TAE. 

As we have seen, there exist a lot of other interesting 
solutions to the TRE, TAE and RWE, listed above, for 



N+l 



A d ■ 4- 
a r ± — — — ia a u> ± 



R%(r)=T> N R%*(r), 



p(r-r+) 



p (r — r_) 



(IV. 3a) 



dr 



R%(r)^(r-r + ) st - 1 ^/ 2 (r-r^y^^/ 2 X N (z ± ),(lV.3c) 



N 



(±l) iV + 1 (-l) 



(JV+l)/2 



(IV. 3d) 
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It is easy to check that in the case of Eqs. re- 
lations pVTTj) and (jTVT3|) produce the standard TSI ifO]) 
and p.6j) with Starobinsky's constants £jv and ®n. 

To reach this result one has to apply the new nontriv- 
ial hidden symmetry of the parameters of (5/v-conflucnt 
Heun's functions, considered as functions of the spin- 
weight s in the special case JV+l = 2|s|. Then: 

{aJ( S ),/3J( S ),7i(s),^( S ),^(s)} = 
= {a ± (-s),/3 ± (- S ), 1± (- S ),S ± (-s),r)±(-s)}. (IV.4) 



Equation (|TO)) follows from (fTTTT]) . (fTTT4|) and (jTT3]) . if 

N+l = 2\s\. As a result, in the case X*(s,z±) = 

Xn(—s, z±) is a solution to the same TRE or TAE but 
with spin- weight (— s). The last fact was discovered for 
the first time in |l|], as a special property of the par- 
ticular solutions to the TRE or TAE, considered there. 
As seen, we obtain a natural explanation of this phe- 
nomenon using the symmetry properties of (5/v-conflucnt 
Heun's functions. 

Now it remains only to point out that 

A) From Eqs. (llV.lj) for cr a = <jp = <r 7 = — a and 
n = s — k one obtains 



k cot ( 



when a = +1: or 



L, S+(9) = +B S(0), S+*(9) =- s S(0) 



d e ,+ -kcot9 = n tf, S+(9) = _,S(0), S+*(9) = +s S(9) 

when a = — 1. Hence, for these particular values of the 
parameters Eq. (|iV.2j) coincides with the standard an- 
gular TSI (HI]). 

B) From Eqs. pV.3[) for <j a =op = a 1 = — a one obtains 

d r .+ = D, R + N (r) = A s +S i?(r), i?+*(r) = _,fl(r), 

and Dn = when a — +1; or 

dr,+ = I»t, R+ (r) = _ a J2(r), i?+*(r) = A s +s R(r), 

and Dn = s£>* when a = —1. Hence, for these particular 
values of the parameters Eq. pV.3a[) coincides with the 
standard radial TSI (|l6j) . 

To our knowledge, the TSI for the other seven cases 
- Eq s. di ll), (IHT31) . (ITlLlal) . (lHI4b1 . (THL511 . (lHI6aTl . 
and pil.6bp have not been studied up to now. In partic- 
ular, relation (|IV.3c[) extends to all cases pil.3[) . (|III.4[) 
the transformation invented for the case 1[) in [l|, [2] . 

3. For the RWE (in units 2M = 1 =>• A = r(r - 1)): 



A 



(N+l)/2 







dr,0,± 

d 



N+l 



R±Jr)=® N>0 R%%(r), (IV.5a) 



Vo,± — -; «cr Q w 

ar 



^(s^T^is-tu;), (IV.Sb) 



Sjv.o = (±l) w+1 $Ptf- 



(IV.5c) 
(IV.5d) 



For the TSI pV3)l we have iV+l = |s|. 

4. The case of the ZE is slightly different. In it we may 
proceed in the following way. 

The solutions Z^(r) to the Zerilli equation are known 
to be simply related with the solutions i?^ (r) to the RWE 
[J] . We write down the corresponding relations (in units 
2M = 1) in the following form (see the appendix): 

iu,R±(r) = e*W± (e-»»4H) , (IV.6a) 

iQZt{r) = ^l-^e-^( r )-^-(e^Mi?±(r)).(IV.6b) 

Let us introduce the notation (r) for Sn solutions 
to the ZE and Z^*(r) for corresponding associate solu- 
tions, and the new operator z rf r> o,± = ^r.o,± + Aq . Then 
after substitution of pV.6a[) into Eq. pV.5a[) we ob- 
tain the following nonstandard Teukolsky-Starobinsky- 
like identities for the solutions to the ZE: 



r+i / , \N+1 / 
2 (*<*r,0,±) f 



d 
dr 



2); 



(IV.7) 



V. A NEW METHOD FOR CALCULATION OF 
STAROBINSKY'S CONSTANTS £ N AND D N 



Equations pV.3d|) . pV.ldj) and pV.5dj) show that in 
all cases the corresponding Starobinsky constants coin- 
cide (up to known numerical factors) with the coefficients 
v n (a, f3,j,S,rj) in the Taylor series expansion p,13j) of 
(5/v-conflucnt Heun's function evaluated for proper values 
of the parameters. These coefficients can be determined 
by the three-terms recurrence relation [l3l . fl5| : 

A n v n = B n v n -i + C n v n - 2 , n=l, ...,oo, (V.l) 

with the initial condition u_i = 0, vq = 1. For 5n~ 
confluent Heun's functions in (|V. 1|) we have coefficients 



K 

Br, 



t 

n 

K + 1 K — /i 



(V.2) 



C n = (n-N-2) 



Relations (|V. 1 1) , (IV. 2|) give a very efficient new method 
for calculation of Starobinsky's constant for all values 
of the spin |s|, confirming the known results and easily 
solving the problem for all cases (|III.1[) - (|III.6|) of validity 
of the (^-condition. 



VI. CONCLUSION 

In the present article, we showed that TSI exist for 
a large class of physical problems when the background 
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vacuum metric is of Petrov type D. We constructed ex- 
plicitly the TSI for different perturbations of any spin 
to the Kerr metric, as well as to the Schwarzschild one. 
Our derivation of TSI is valid for complex values of the 
separation constants w and E, as well as for complex val- 
ues of all other parameters. The obtained results can be 
applied to both the black holes and naked singularities. 

The exact solutions to the TRE, TAE, RWE and ZE 
were briefly described in terms of confluent Heun's func- 
tions. One can find more details in the articles [T2|. Hav- 
ing in hands these explicit solutions we are able 

1. To write down the corresponding general solutions. 

2. To formulate different boundary problems. 

3. To derive all properties of the solutions using ex- 
plicitly the properties of Heun's functions. 

In particular this article accomplished the following 
aims: 

i) We found the TSI and all their possible general- 
izations have a common origin - the recently unveiled 
property (|I.12[) of the confluent Heun function, related 
with the 5n condition (|I.9|) . see [l3T| . As a result of this 
condition, the multiple derivative of order (N + 1) of a 
given (5/v-confluent Heun's function becomes proportional 
to another (associated) confluent Heun's function. The 
coefficient is precisely (7V+1)! times the (7V + l)-th coeffi- 
cient in the Taylor series expansion of the first function. 
This algebraic property is universal and not related to 
the usually considered specific boundary problems, nor 
with other properties of their solutions like regularity, 
integrability, etc. Using the relation of the exact solu- 
tions at hand with the confluent Heun's function, one 
can express this property directly in terms of the very 
solutions. Thus one obtains the TSI for the TRE, TAE, 
RWE and ZE in a uniform and universal way. 

ii) We discussed all classes of S n solutions to the above 
problems. For the specific class of (5Ar-solutions, defined 
by relation 1|) . the novel symmetry pV.4[) of the pa- 
rameters of <5Ar-confluent Heun's functions jL8| ensures 
that its multiple derivative of order (N + 1) is a solution 
to the same Heun equation. Thus, this symmetry yields 
a second solution to the same TME equation. As a re- 
sult, we obtain the known pair of solutions to the TME 
with the opposite signs of their spin- weights (±s) 

iii) We extended the familiar results for TSI to all 
possible cases on a firm and general mathematical ba- 
sis. Thus, using the <5/v condition we found the possible 
generalizations of the TSI. 

iv) We made transparent the origin of some of the prop- 
erties of the solutions to the TRE, TAE, RWE and ZE, 
relating them with the specific form of the parameters in 
confluent Heun's functions, which solves these equations. 
For example, such unexpected property is the symme- 
try pV.4p , which is specific for large class of solutions to 
TRE and TAE. This symmetry does not take place for 
all confluent (5/v-Heun's functions. 

v) Using the recurrence relation for the Taylor series 
expansion of the confluent Heun's function we proposed 
a new efficient and universal method for the calculation 



of Starobinsky's constants. 

Thus, the confluent Heun function was proven to be an 
adequate tool for solving of the above problems. The us- 
age of this function gives a natural treatment and deeper 
understanding of the Teukolsky-Starobinsky identities, as 
well as of the solutions to the Teukolsky equations. The 
results, presented here, are a basis for future physical 
applications. 
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APPENDIX A: ON THE RELATION BETWEEN 
THE SOLUTIONS TO THE RWE AND ZE 

The relation between the solutions to the RWE and 
ZE were discovered by Chandrasekhar, see for example 
[3] and the references therein. We slightly modify the 
form of these relations using the notation Z^(r) for so- 
lutions to the Zcrilli equation and Rq(t) for solutions 
to the RWE. These functions have additional constant 
phases in comparison with the functions Z^\r), used in 
[J] (Note that the meaning of the signs '±' in our nota- 
tion is completely different see [Uj].) Our choice of these 
additional phases ensures the appearance of the common 
factor iuj simultaneously in the two of the Eqs. (|IV.6I) . 

As usual, here we consider only the ZE for gravitational 
waves with spin |s| = 2=>A+l = 2. Then in units 
2M = 1 we obtain 

e * w = k-ir^(i + (| _^ +2)r )"' (A.i) 

ui = l(l-l)l(l+l)(l + 2), 1 = 2,3,..., 



U! = y uj 2 + ujf ^ for uj ^ ±iu>i . 

The Eqs. PV.6|) do not relate the solutions Z Q (r) and 
i?J(r) when uj = ±iuj[. Instead, in this case one obtains 
from the Eqs. pV.6[) the well known algebraically special 
solutions in the following explicit quasipolynomial form: 



7 



R$(r) = const x e +^ r \ (A.2a) 
Z±(r) = const xe^ r ). (A.2b) 

Thus we see that the case of the algebraically special 
solutions (|A.2|) is a degenerate one. 

The system of the ordinary differential equations 
(|IV.6j) can be easily split to the decoupled RWE and ZE, 
written in the form 

R" + uj 2 R = (-ip" + V>' 2 ) R, (A.3a) 
Z" + Cj 2 Z = (+^j" + ^Z. (A.3b) 



Here the prime denotes differentiation with respect to 
the tortoise coordinate r* = r + In |r — 1|. This splitting 
reflects the physical independence of the axial and polar 
perturbations to the Schwarzschild metric. 
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